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0.0 ABBREVIATIONS 


CORALA 

G.S.A. 

IR 

M/S 

RH9100 

S.S.D.A. 

VR 


Correlation iVnalysla 
Grey Scale Adjustment 
Infrared .85 micron 
Multlapectral ratio 
RAKTFK 9100 

Seqviential Similarity Detection Algorithm 
Visual Red .65 micron 
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I . FOREWORD 


This volume of the "Onboard Utilization of Ground Control Points 
for Image Correction" final report provides a detailed description of 
the software hat was utilised for image correction accuracy measure- 
ment. This soitware tvas developed under an independent research task 
DllR "Video Guidance Technology" and is included here for complocencss . 

Three other volumes have been Incorporated into the final report. 
Volume 1 provides an executive summary, Volume II contains a detailed 
description of the study and simulation results, and Volume ill is 
an appendix describing the software developed for simulation of the 
onboard navigation system. 
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1 . 0 INT RODU CTION 

The correlation analysis program Is written to allow the user various 
tools to anuly/o different ccrrelatlon algorithms. These algorithms were 
tested using L.\NDSAT Imagery in two different spectral bands. Image data 
of the same area was taken at two different times. In the Imagery taken a: 
tine t, a reference landmark la extracted by the user. An attempt Is then 
made to find the corresponding landmark in the image data taken at time t 2 
by using the various correlation algorithms. 

Each correlation performed returns various monitoring values. The 
time required for the correlation is returned as well as the correlation co 
efficient. Cloud percentages in both landmark and search areas are dis- 
played. The absolute beat fit element and line values are also calculated. 

The user controls the flow of the CORALA programs by defining landmark 
sizes, search area sizes, algorithm types, and which bands to use for the 
correlation run. Currently, the maximum search area size Is 6A x 6^ pixels 
and the maximum landmark size is 56 56 pixels. Correlation may be per- 

formed on the imagery in any of three display areas. The screen setup on 
the RM9100 is lllustrateu in Figure 1.0. In this example, the user has 
defined a landmark in band 1. The correlation program would attempt to 
find this l^'.dmark In band 1 of the data taken at time t 2 * 

Currently, three classification algorithms have been implemented. 

Table 1.0 describes two of the three algorithms. The third algorithm, 
S.S.D.A. with thresholding, is a modification of the S.S.D.A. algorithm. 
This algorithm will proceed as long as the current correlation value being 
computed Is less than the previous minimum. 
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The primary control varlablca which can be altered for each particular 
run by the uaer Include; 

- Scene location 

- GCP location 

- Search area size 

- GCP size 

The observables on which the analysis is performed Imlude: 

- Time of correlation 

- Rate of convergence 

- Percentage difference in the rate of convergence for live-scene 

correlation versus autocorrelation 

- Probability of false lock 

- False lock deevetion feasibility 

- Ease of Implementation In a hardwired system 

- Location of correlator lock 

The quantitative results of the analysis are provided In a hard-ropy print- 
out, and a correlation surface plot Is displayed on the video monitor 
(Figure 1.0,1). These data provide the basl« of a tradeoff between the 
various algorithms. 


Li 
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2 • 0 USE R’S GUIDE 
2 • 1 Introdu ction 

The Correlation Analysis Program, CORALA, has two defined user inodes; 
they are the Interactive and the batch modes. The Interactive mode queries 
the user for Information concerning the current running of the program. 

Using the batch mode, this data has already been prestored, and all the user 
must do Is Initiate the running of the C program. 

An auxiliary task called ALTOCOR, aut.jmatlc correlation, Is used by 
the user when batch mode Is the desired mode of operation. This task basi- 
cally duplicates the query sect:* on that would be required for an Interactive 
CORALA run. 

Under either mode of operation, the results of the correlation may be 
stored in a data base file containing previously-stored runs. The informa- 
tion contained in the data base is described in Section A. When the data 
base storing option Is selected, a predefined run sequence of landmarks, 
sizes, and algorithms are used. 

The following types of runs may be performed: 

(1) INTERACTIVE MODE, NO DATA BASE STORAGE 

(2) INTERACTIVE MODE. DATA BASE STORAGE 

(3) BATCH MODE, DATA BASE STORAGE 

CORALA is initiated by the user by the following command sequence: 

Log in to tlie PDP 11/70 
ASN DRl: ■ SY : 

SET/UI^ - ^[ll2,5''| 

MOU DR l : LAN DSAT 
PIP * .DAT /UN/NM 
RUN CORALA 
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From this point on. th< us«r Is promptsd dcpsnd^^nc upon cha cyps of 


run being performed. 

Double underline indicates user input. 

Types of Runs 

2.2.1 T YPE 1 - INTERACTIVE MODE. NO DATA BASE STORAGE 

The following is the sequence of questions asked the user when 
under a Type 1 run; 

(1) Will this run use a prestored run controller (Y or N) N 

Inputting an N will cause this to be in the interactive 
mode . 

(2) Is this run to update the corr. data base (Y or N) N 

This will put the user in the no data base storage mode. 

(3) Input tet(t case number; nvmibe£ 

This number may be ±32767. 

(4) Input file name for search area 1; name B2.ext. 

Each image file name is designated by a name appended 
with a B2 or B4 signifying Band 2, .65 micron, or Vand 4, 
.85 micron. For this question, B2 is appended. The file 
extension, ext, is a 3-letter file descriptor type, 
usually DAT. 

(5) Input starting element, line, and ratio (315): n,n ,n 

Each of there numbers may be in the range ±32767. De- 
fault display ratio is 1. The element and line inputs 
define the absolute location of the imagery to be 
displayed. 

(6) The user Is prompted to flip switch 1 of the Joystick up. 

Then the user positions the cursor, by using the Joystick, 




to the center of Che desired search area. Once pcsi cloned, the 
user must flip switch 4 up, Chen down. 

(7) Enter option for third display area — 

1 ■ MulclspecCral ratio 

2 ■■ Grey scale adjustment 

Option type: 1 or 2 

This defines the type of enhancement used for the third 
display area. 

If a 2 has been input for question 7, the following question 
ic asked: 

(7a) Input the type of G.S.A. to use — 

1 « Standard deviation method 

2 “ Histogram equilization method 
OPTION >> I or 2 

This defines the type of grey scale adjusCment needed. 

(8) Enter element, line offset between the 2 setnes (214): n . n 

Each of these n'-irabers may be in the range of ±32767 . 

(5) Input file name for search area 2: name B2.ext 

See Instruction (4) . 

(10) Input algorithm number(8) to use 
SSDA 1 

CLASSICAL 2 

SSDA with threshold 3 

Algorithm(s) type: n l,n,..| 

The user may input one, two, or all three correlation 
algorithms to perform the correlation. 


I 
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(11) Input display aroa(s) covrolatlon sliould take place 

1 • .65 micron band 

2 - .65 micron band 

3 “ Multlspectral or C.S.A. area 

Display area(s) : n | , n . . . [ 

The u8f«r may perform the correlation in any or all bands. 

(12) Inp’Jt the landmark sizes to be used. The largest size must 
appear first In the list. Maximum l.ra. size - 54 

Landmark sizes ; n |,n...| 

Up to ten .'andmark sizes may be input. 

(13) Input 1 .m. sizes to nub-plxel registered: n j,n...] 

Tlio l.m. sizes in question 12 can be used. 

(14) Input l.m. sizes to have their corr. surfaces plotted: n|,:i...| 

The l.m. sizes in question 12 can be used. 

(15) 3'he user is then prompted to locate the landmark center using the 
same tn.struct ions as question 6. 

From here on. all correlation is autonomous. The us. r can ait back. 
Wlicn the program finishes, printout of the run may be received by 
tvping in - IMP COILVLA.DAT/SP 
TYFF 2 - INTERACTIVE MODE, DATA BASE STORAGE 

The following Is the sequence of questions asked the user when 
under a tvpe 2 run: 

(1) Will this run use a prestored run controller (Y or N) N 

Inputting an N will cause this to be in the interactive 
mode . 

(2) Is this run to update the corr. data base (Y or N) Y 

This will put t’,\e user in the data base storage mode. 





(3) Input file namp for search area 1 ; name B2.ext 

Each image tile name ia designated by a name appended with 
a B2 or B4 signifying Band 2, .65 micron, or Band 4, .85 
micron. For this question, B2 is appended. The file ex- 
tension, ext, ia a 3-letter file descrlptir.i type, ULi»'.ally 
DAT. 

(4) Input starting element, line, and ratio (315): 

Each of these nunbers may be in the range ±32767. Default 
display ratio is 1. The element and line inputs define 
the absolute location of the Imagery to be displayed. 

(5) The user is prompted to flip switch 1 of the joystick up. Then 
tlie user positions the cursor, by using the joystick, to the 
center of the desired search area. Once positioned, the user 
must flip switch 4 up, then down. 

(6) Enter element, line offset between the 2 scenes (214) : n ,n 

Each rf these numbers may be in the range of ±32767. 

(7) Input file name for search area 2 : name B2.ext 

See Instruction (4) . 

From here on, all correlation is autonomous. The use can sit back. 
Wlien the program finishes, printout of the run may be received by 
typing in > PIP CORALA. DAT/SP 

Default values used for a type 2 run: 

(1) Band 3 uses multispectral ratlolng. 

(2) All 3 algorithms are used for the correlation. 

(3) Bands 1 and 2 are used for the correlation. 

(4') Landmark sizes of 48, 32, 24, 16, and 8 are used. 

(5) Subpixel registration is performed on all l.m. sizes. 

(6) No correlation surfaces are displayed. 
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2.2.3 TYPE 3 - BATCH MODE. DATA BASE STORAGE 

The tollowlng queatloni are aeked Che user when under a type 3 

run; 

(1) Will this run use a presCored run controller (Y or N) Y 

This will cause batch mode to become active. 

(2) Is this run to update the corr. data base (Y or N) Y 

This activates Che data base storage mode. 

From this point on, everythin', Is autonomous. Defaults, landmark 
sizes, etc., are the same as those for a type 2 run. 
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3.0 PROGRAMMER'S SECTION 


3. 1 Introduction 

The CORALi\ programs use a tree structure to define their operation. The 
root of the tree Is the main program CCRALA. This main program acta as .a 
driver to the three major factors: CORINI, CORMID, and CORFIN. These pro- 

grams In turn act as drivers to other subroutines. The tree structure may be 
seen In Section 3.4 (Figure 3.4b). 

All programs were written in FORTRAN 4 PLUS with the exception of the 
I/O routine /VDJUST. These programs were designed and Implemented on a DEC 
PDP 11/70 running RSX ll-M operating system. The graphics device used was a 
RAOTEK 9100 320 x 256 pixel resolution graphic display system. 

The C0R/U.A subroutines are classified into three types: 

O) Progriira specific 

(2) Library routines 

(3) I/O routines 

3 . 2 Variable Conventions 

Wlienever possible, the following variable naming conventions were used: 

(1) Any variable beginning with an I Is considered an 1*2 variable (2 bytes) 

(2) Any variable beginning with an R is considered an R*4 variable (4 bytes) 

(3) Anv variable beginning with an A or L is considered byte variable 

(1 bvte) . 

(43 Anv variable not following the previous conventions are declared using 
standard F0RTR(\N variable conventions. 

3 . 3 Common Blocks 

Common blocks are used extensively in the C0RAL<\ software package. 

These common blocks facilitate the sharing of common variables. The blocks 
of code are inserted in the appropriate source programs by using the INCLUDE 
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statement. This eliminates the possibility of a typographical error when 
typing in multiple common blocks. 

3.3 1 DATA.COM 

This common block contains all scratch arrays and various other 
data arrays and variables. This is sort of a catch-all common block. 
COMMON /DATA/ IWIN(64,64) ,ILM(64 ,64) ,IBUF(320) ,ISASZE,IJ0Y(4) , 

IESL(3) ,ILSL(2) ,IBFE,IBFX,ITPFE,ITBFT,ICLDTH, 

LBUF( 5 12), IRATIO , IBOPT , IGSATY , RTHRES , VEGBAR , BARWAT , 
RBFE,RBFL,IRED(256) ,IRED1(256) ,IC0NT 
INTEGER*: TWIN , ILM , IBUF , ISASZE , IJOY , lESL , ILSL , IBFE , IBFL , 

TTBFE , ITBFL , ICLDTH , IRATIO , IBOPT , ir^ATY , IRED , 

ICONT.IREDl 

REAL* 4 RTHRES, RBFE.RBFL, VEGBAR, BARWAT 

LOGICAL*! LBUF,LUBF1(8192) 

EQUIVALENCE (LBUFl ( 1) , ILM( 1) ) 

IWIN — SCRATCH ARRAY 

ILM — SCR.'vTCH ARRAY 

IBUF — SCRATCH ARRAY 

ISASZE — SEARCH AREA SIZE IN PIXELS 

IJOY ~ JOYSTICK SCRATCH BUFFER 

lESL — ELEMENT START ARRAY 

ILSL — LINE START ARRAY 

IBFE — BEST FIT ELEMENT 

IBFL — BEST FIT LINE 

ITBFE — TEMPORARY BEST FIT ELEMENT 

ITBFL — TEMPORARY BEST FIT LINE 

ICLDTH — VR. CLOUD THRESHOLD VALUE 
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LBU? — SCRATCH ARRAY 
IRATIO — DISPLAY RATIO 

IBOPT — OPTION TYPE FOR THIRD DISPLAY AREA 
IGSATY — OPTION TYPE FOR GREY SCALE ADJUSTMENT’ 

RTHRES — CLOUD PERCEi.TAGE MAXIMUM VALUE 
VECBAR — VEGETATION/BARE EARTH BOUNDARY RATIO 
BARWAT -- BARE EARTH /WATER BOUNDARY RATIO 
RBFE — BEST FIT ELEMENT (REAL) 

RBFL — BEST FIT LINE (REAL) 

I RED — SCR/.TCH ARRAY 
IREDl — SCRATCH ARRAY 
ICONT — AUTO CORREL^\TION FLAG 
3,3.2 DATABS.COM 

This common block contains the variables and arrays for the data 
base storage mode. 

COMMON /DATABS/ lOUT ( 1 792) , IAUT0( 81) , IACT( 8 1 ) , IDBUP , IRECNO 
LOGICAL*! LOUT(3584) 

REAL*4 ROUT (896) 

INTEGER*2 lOUT , lAUTO , lACT , IDBUP , IRECNO 

EQUIVALENCE (lOUT(l) .LOUT ( 1 ) . ( lOUT ( 1 ) ,ROUT(l))) 

lOUT — OUTPUT ARRAY USED TO HOLD ALL CORRELATION PERTINENT DATA 
I AUTO — ARRAY CONTAINING THE AUTO CORRELATION VALUES 

lACT — ARRAY CONTAINING THE ACTUAL CORRELATION VALUES 

IDBUP — DATA BASE UPDATE FLAG. 'Y '-UPDATE, OTHERWISE NO UPDATE 

IRECNO — CURRENT NUMBER OF TEST CASES STORED IN THE DATA BASE 
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3.3.3 FILES.COM 


This common block contains the LANDSAT file names end the sizes 
of these names of the Images curre.itly displayed. 

COMMON /FILES/ IF2SZE,FNAMEl (20) ,FNAME2(20) 

LOGICAL*! FNAMEI .FNAME2 
INTECER*2 IF1S2E,IF2SZE 

IFISZE — SIZE OF THE NAME FOR THE FIRST DATA FILE 
IF2SZE — SIZE OF THE NAME FOR THE SECOND DATA FILE 
FNAMEI — ARRAY CONTAINING THE FIRST FILE NAME 
FNAME2 -- ARPJiiY CONTAINING THE SECOND FILE NAME 
3.3.4 HELiVDER.COM 

This common block contains all the HEADER arrays used for the 
extraction and display of the LANDSAT image HEADER block. 

COMMON /HEADER/ SID( II) ,E.XP0DT(7) ,F0RCEN(14) ,SUNANG(21) ,SUNANA(3) , 
NADIR(14) ,HEADNG{3) 

LOGICAL* I SID, EXPODT , FORCEN , SUNANG , SUNANA , NADIR , HEADNG 

SID — SATELLITE IDENTIFICATION 

EXPODT — EXPOSURE DATE 

FORCEN — FORMAi CENTER 

SUNANG — SUN ANGLE 

SUNANA -- SUN AZIMUTH ANGLE 

NADIR — NADIR LATITUDE AND LONGITUDE 

HEADNG — SATELLITE HEADING 
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3.3.5 IO.COM 


Thla common block concalna all pertinent I/O variables used to 
talk to the RM9100. 

COMMON /lO/ I0ST(4) .IPARAM(6) ,IRMCMD(30) ,IDS 

lOST — I/O STATUS ARRAY 
IPARAM — INSTRUCTION PARAMETER ARRAY 
IRMCMD — RAMTEK INSTRUCTION ARRAY 
IDS — DIRECTIVE STATUS 
3.3.6 OPTS ■ COM 

This common block contains the user Input control options. 
COMMON /OPTS/ ISURF(IO) .IBOX(IC) ,IALG(3) ,IBANDS(3) .ISBOX(IO) , 

INBOX, INBAND, INALG 

INTEGER* 2 ISURF , IBOX , I ALG , I BANDS , I SBOX , INBOX , INBAND , INaLG 

ISURF — CONTAINS THE LANDMARK SIZES TO HAVE THEIR CORRELATION 
SURFACES DISPLAYED 

IBOX — CONTAINS THE LANDMARK SIZES TO USE 
lALC — CONTAINS WHICH ALGORITHMS TO USE 
IBANDS — CONTAINS WHICH BANDS TO USE 

ISBOX — CONTAINS WHICH L.M. SIZES TO SUBPIXEL PJIGISTER 
INBOX — THE NUMBER OF LANDMARK SIZES INPUT 
INBAND — THE NUMBER OF BANDS INPUT 
INALG — THE NUMBER OF ALGORITHMS INPUT 
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3.3.7 COMMON BLOCK REFERENCE LIST 



DATA 

DATABS 

FILES 

^HEADER 

10 

OPTS 

ADJ 

X 






AUTO 

X 






AUTOCR 

X 

X 





aouDc 

X 






CLOUDS 

X 






CLSCOR 

X 






CORAI.A 

X 

X 

X 

X 

X 

X 

CORFIN 

X 

X 




X 

CORG 

X 






CORINI 

X 

X 

X 

X 



CORMID 

X 

X 




X 

CORSDA 

X 






CRALGS 

X 

X 





CRCORR 


X 





DISPLE 

X 

X 





DIVIO 

X 






GETHUR 




X 



GRID 

X 






GSAHIS 

X 






GSASTD 

X 






LMDISP 




X 



NXDREC 

X 






POS 

X 






RESET 





X 


RMCPCK 





X 


RMERSE 

X 






RMRDIM 





X 


RlfTEJCT 





X 


RMWRIM 





X 


SAREA 

X 

X 





SCALE 

X 






SPIX 

X 






WIO 





X 


WRTTSC 

X 

X 

X 
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3.4 CORALA Program Tr'ge Structure 


CORALA 



Figure 3.4a 
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3.5 Program upaclflc routinaa art thoaa rouclnaa uaad apaclflcally for cor- 
ralatlon analysla. Thaaa routinaa are. 




SECTION 

% 

i • 

MXJ 

3.5.5 

2. 

AUTO 

3.5.6 

3. 

AUTOCR 

3.5.7 

4. 

CLOUDC 

3.5.8 

5. 

CLOUDS 

3.5.9 

6. 

CLSCOR 

3.5.10 

7. 

CORALA 

3.5.1 

8. 

CORFTN 

3.5.4 

9. 

CORf 

3.5.11 

10. 

COR IN I 

3.5.2 

11. 

CORMID 

3.5.3 

12. 

CORSDA 

3.5.12 

13. 

CORSMN 

3.5.13 

14. 

CRALGS 

3.5.14 

15. 

CRCORR 

3.5.15 

16. 

DISPLE 

3.5.16 

17. 

DISPLY 

3.5.17 

18. 

DIVID 

3.5.18 

19. 

GRID 

3.5.19 

?,0. 

GSAHIS 

3.5.20 

21. 

GSASTD 

3.5.21 

22. 

SAREA 

3.5.22 

23. 

SCALE 

3.5.23 

24. 

SPIX 

3.5.24 

25. 

SPLPRO 

3.5.25 

26. 

WEIGHT 

3.5.26 
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3.5.1 CORALA (Correlation Analysis) 

CORALA is the main program for the correlation analysis pro- 
grams. This program will allow the user to run the correlation 
programs either in a semi-batch mode or an interactive mode. 

CORALA assigns the Joystick to unit 2, resets the RM9100 and 
opens the data base and output files. This program then drives the 

major subroutines CORINI, CORMID, and CORITN. 

CORALA 


RESET THE RM9100 


SET THE JOYSTICK TO UNIT 2 


OPEN THE OUTPOT DATA FILE 




IS THIS A BATCH RUN? 


OPEN THE BATCH 
CONTROLLER FILE 


NULL 




IS THIS A DATA BASE UPDATE RUN/ 


GET THE CURRENT 
RECORD NUMBER 


NULL 

CALL CORINI 


CALL CORMID 


CALL CORFIN 


DELETE THE CLOUD TEMPORARY FILE 

IS THIS A DATA BASE 

UPDATE 

UPDATE THE CURRENT 
NUMBER OF TEST RUNS 


NULL 

IN THE HEADER 




EXIT IF THIS IS NOT A BATCH RUN OR 
IF THE BATCH CONTROLLER IS OUT OF DATA 


CLOSE ALL FILES 


RESET THE RM9100 
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3.3.2 CORINl (Correlation Initialization) 


The CORINl r-^utine is the initialization subroutine tor the 
correlation analysis programs. This subroutine will display the 
image data in display arean 1, 2, and 3. These d' lay areas will 
contain the usei -defined search area tai'en at two different times. 


RM9100 


.65p 


.85p 


M/S or 
G.S.A. 




. 65u 


.85p 


M/S or 

G.S.A. 





data taken 
at time Tl 


data taken 
at time T2 


This routine will call otlier routines that display two 320 x 
120 pl.xel image areas of the same imagery in two different bands. 
The use then defines the search area size and defines the center of 
this search area in the imagery displayed. The user will then de- 
fine which type of algorithm for display area 3 should be used. 

Tlie corresponding data taken at time 2 will then be displayed. 
Control is Chen returned to CORALA. 
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CORINI 


IS THIS A DATA BASE RUN? 


GET THE TEST CASE NUMBER 
FROM THE DATA BASE 


GET THE TEST CASE NUMBER 
FROM THE USER 


IS THIS A BATCH RUN? 


GET THE FILE NAME FOR 
THE TIME I DATA FROM 
THE BATCH CONTROLLER 



GET THE FILE NAME FOR 
THE TIME 1 DATA FROM 
THE USER 


IS THIS A BATCH RUN? 




GET THE STARTING ELEMENT, 
LINE, RATIO FROM THE 
BATCH CONTROLLER 


GET THE STARTING ELEMENT,! 
LINE. RATIO FROM THE 
USER 


DISPLAY THE 2 SPECTRAL BANDS OF DATA 


GET THE SEARCH AREA MD SIZE FROM THE USER 


GET THE THIRD AREA DISPLAY TYPE AND DISPLAY 


7! 

!!• 


IS THIS A BATCH RUN? 


GET THE SECOND FILE NAME 
FROM THE BATCH CONTROLLER 

GET THE SECOND FILE NAME 
FROM THE USER 

DISPLAY THE LOWER 
WITH TIME i 

H/'iLF OF THE SCREEN 
IMvXGE DATA 

CREATE THE CLOUD DATA FILE SPECIFYING WHICH 
PIXELS IN THE SECOND SET OF IMAGERY ARE CLOUDS 
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3.5.3 CORMID (Correlation Middle Routine) 


The CORMID subroutine queries the user to input certain control 
parameters. In a data base storage mode, these parameters are pre- 
defined. The questions asked are: 

(1) Algorithm(s) to be used for correlation 

(2) Band(s) for the correlation to be tested 

(3) Landmark size(s) to be tested 

(4) Landmark size(s) to be subpixel registered 

(5) Landmark size(s) to have correlation surface plotted 

The default values used for a data base storage mode are: 

(1) S.S.D.A. and classical correlation 

(2) Bands 1 and 2 

(3) Landmark sizes of 48, 32, 24, 16, and 8 

(4) Subpixel landmark sizes of 48, 32, 24, 16, and 8 

(5) No correlation surfaces plotted 


FORM/\T : C^\LL CORMID 


CORMID 


TNIT1;\LIZE COUNTERS 


IS THIS A DATA BASE UPDATE? 


SET UP DEFAULT 
VALUES 


QUERY USER FOR 
VALUES TO 
BE USED 
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3.5.4 CORFIN (Correlation Flniah) 


Once all the Input values and Image data are received from the 
user, CORFIN proceeds to perform the actual correlation analysis. 
This subroutine drives the analysis programs. 

Hard-copy correlation results are printed from this routine. 
This routine also stores the correlation results Into the datn base 
If this Is a data base mode run. 


FORMAT ; C,\LL CORFIN 
CORFIN 


IS THIS A BATCH MODE RUN? 

GET l.v:d>l\rk loc.\tijn 
FROM THE BATCH FILE 
1 . 1 

GET LANDMARK LOCATION BY ALLOWING 
THE USER TO USE THE JOYSTICK 


PERFORM AUTO CORRELATION 


DRAW A BOX AROUND THE REFERENCE L/INDMARK 


SET UP TIMER TO ZERO AND START TIMING 


CALL THE ROUTINE TO PERFORM CORRELATION 


STOP TIMER AND STORE THE VALUE 


DETERMINE THE L.VNDMARK CLOIFD PERCENTAGES 


DETERMINE THE CORRELATION COEFFICIENT 


DISPLAY THE CORRELATION SLTIFACE IF DESIRED 


PERFORM SUBPIXEL REGISTRATION 


RESTORE SCREEN 


OUTPUT DATA 


DO FOR THE LANDMARK SIZES SPECIFIED 


DO FOR THE NUMBER OF ALGORITHMS SPECIFIED 


DO FOR THE NUMBER OF B.\NDS SPECIFIED 
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3.5.5 ADJ (Adjust data) 

ADJ is the driver routine used for grey scale adjustment of 
image data. Two methods of grey scale adjustment are currently 
implemented. One method is the standard deviation method. The 
other is the histogram equilizatlon method. ADJ allows the user to 
specify the appropriate method and then calls the subroutine that 
implements this method. 

Two flags are passed to this subroutine. IFLAG is a decision 
flag declaring which method of grey scale adjustment is to be used. 
If IFLAG equals zero, the user is prompted for the G.S.A. type. 
IFLAG2 is another decision flag used for hard-copy output. This 
variable will be used by the adjustment routines. ILS is the line 
start of the imagery. 

FORMAT: CALL ADJ (ILS , IFLAG, IFLAG2) 

ADJ 


IS IFLAG i 0 


NULL 



HAVE USER INPUT 
IFLAG 

IS IFLAG - 2 

CALL THE ROUTINE TO 
PERFORM HISTOGRAM 
EQUILIZATION G.S.A. 

CALL THE ROUTINE TO 
PERFORM STANDARD 
DEVIATION G.S.A. 
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3.5.6 AUTO (Auto Correlation) 

AUTO la used In conjunction with the subroutine AUTOCR to per- 
form nuto correlation. This subroutine will be called 81 times by 
AUTOCR to create a 9 x 9 element auto correlation surface. 

This routine performs the correlation algorithm on one place- 
ment. The algorithm type is defined by ITALG. ITAUTO Is the 
correlation value determined by this placement. IX and lY are off- 
sets In the array TWIN which contains the landmark Image data. 
RLMEAN is the landmark mean ani. ISIZE is the landmark size. ICORR 
is used by the S.S.D.A. with thresholding algorithm and contains 
the current minimum correlation value. 

FORMAT: CALL AUTO ( ITAUTO , IX, lY , RLMEAN , ISIZE , ITALG .ICORR) 

AUTO 
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3.5.7 AUTOCR (Auto Correlation Driver) 


AUTOCR Is the auto correlation driver routine. Auto correla- 
tion Is the process of correlating a defined landmark upon Itself 
using a ^-element offset In all directions. What this will do Is 
describe the best possible correlation surface for this landmark. 

The placement values ’•e stored In a 9 x 9 element array called 
I AUTO. 

ISTE and ISTL are the starting screen coordinates of the land- 
mark. ILMSZE Is the landmark size and ITALG is the classification 
algorithm being used. 

This subroutine will store the auto correlation values into the 
data base If the data base option is defined. 

FORMAT: CALL AUTOCR( ISTE, ISTL, ILMSZE, ITALG) 

AUTOCR 


READ IN THE LANDMARK DATA AS WELL 
AS A 4-PIXEL BUFFER ON ALL SIDES 


DETERMINE THE MEAN 


DETERMINE THE AUTO CORRELATION 
VALUES BY CALLING AUTO 


IS THE DATA BASE UPDATE FLAG ON? 


STORE IN THE DATA BASE 
THE AUTO CORK. VALUES 


NULL 
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3.5.8 CLOUDC (Create Cloud File) 


This subroutine takes a plxel-by-plxel ratio of two bands of 
one scene and creates a cloud classification array and stores this 
array in the data file CLOUDS.TMP. If a pixel is classified as a 
cloud, the sign bit is set in the cloud array. If this isn't a 
cloud pixel, the array element is set to zero. 


FORMAT : CALL CLOUDC 


CLOUDC 


/ 

READ A LINE OF SEARCH 
AREA DATA 



IS THIS A CLOUD PIXEL? 



SET THE 

SET THE 



CLOUD BIT ON 

BIT OFF 



DO FOR THE SEARCH AREA SIZE 


SAVE THIS LINE IN THE DATA FILE 

DO FOR THE SEARCH AREA 

SIZE 
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3.5.9 CLODPS (Set Cloud Bit On) 

This subroutine will read the cloud claasification data created 
by CLOUDC. The displayed search area data ia contained in IWIN. 

This subroutine will turn on any cloud pixels in this search area 
array. By turning the sign bit on, none of the image data is lost. 
This subroutine is used in conjunction with the correlation algor- 
ithms. This allows for an easy means to distinguish clouds and to 
not use them in the correlation. 

FORMAT: CALL CLOUDS ( IS AVE) 

CLOUDS 


OPEN THE CLOUD DATA FILE 

IS ISAVE - 1 

DELETE THE 
CLOUD FILE 

READ THE CLOUD FILE 
AND TURN ON ALL 
CLOUD PIXELS 

SAVE THE CLOUD FILE 
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3.5.10 CLSCOR (Clasilcal Correlator) 


Thia subroutine calculatea the claaalcal correlator values for 
all landmark placements. The array IWIN contains the search area 
data. The array ILM contains the landmark data. RLMEAN is the land 
mark mean and ILMSZE is the landmark size. The classical correlator 
algorithm is: 


//lines //elements 

I 

i-1 


^ (IWIN/J.I) ~ iwmean)*(ilm(j,i) - ilmean) 


j-1 


//lines //elements 


//lines //elements 


1^1 (iWIN(J.I) - IWMEAN)2*^J^ (ILM(J,I) - ILMEAN)2y 


This subroutine performs an exhaustive search of the landmark 
data versus the search area data. The best fit is where the correla- 
tion value is at a maximum. The generated correlation values have 
their range between +1 and -1. These values are multiplied by 32767 
and stored back into IWIN. 


F0R:iATt CALL CLSCOR (RLMEAN , ILMSZE) 


CLSCOR 


SET THE BEST FIT ELEMENT AND LINE TO I 


DETERMINE THE fF^CH AREA MEAN 


SET THE CLOUD COUNTER AND 
CORRELATION VALUES TO 0 


IS THIS A CLOUD PIXEL? 


INCREMENT 
CLOUD COUNTER 


ADD INTO THE CORRELATION 
VALUE FOR THIS PLACEMENT 


DO FOR THE NUMBER OF ELEMENTS 
IN THE LANDMARK 


DO FOR THE NUMBER OF LINES 
IN THE LANDMARK 


ADD IN CLOUD PIXELS USING CLOUD ALGORITHM 



IS THIS A 

BEST FIT? 

SAVE BEST FIT 

NULL 

' STORE THE CORRELATION VALUE 


DO FOR THE NUMBER OF ELEMENT PLACEMENTS 


DO FOR THE NUi«IBER OF LINE PLACEMENTS 


B-33 
















3.5.11 CORC (CorrelAtlon Using SSDA) 

This subroutlns calculates the S.S.D.A. correlation values for 
all landmark placementa. The array IWIN contains the search area 
data. The array ILM contains the landmark data. RLMEAN ia the 
landmark mean and ILMSZE is the landmark size. The S.S.D.A. 
algorithm la defined: 

''/lines //elements , i 

ill jll *b«UlWIN(J ,1) - IWMEAN) - (ILM(J.I) - ILMEANH 

This subroutine performs an exhaustive search of the landmark 
data versus the search area data. The best fit Is where the correla- 
tion value determined by the S.S.D.A. algorithm Is at a minimum. 

The generated correlation values a."? restored into IWIN. 

FORMi\T: CALL CORG ( RK-IE AN , ILMSZE) 
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CORC 


SET THE BEST FIT ELEMENT AND LINE TO 1 


DETERMINE THE SEARCH AREA MEAN 


SET THE CLOUD COUNTER AND 
CORRELATION VALUE TO 0 


IS THIS A CLOUD PIXEL? 

INCREMENT 
CLOUD COUNTER 

ADD INTO CORRELATION 
VALUE FOR THIS PLACEMENT 


DO FOR THE NUMBER OF ELEMENTS 
IN THE LANDMARK 


DO FOR THE NUMBER OF LINES 
IN THE LANDMARK 


ADD IN 

CLOUD PIXELS US INC CLOUD ALGORITHM 


\ 

IS THIS A BEST FIT? 

X 

SAVE 

BEST FIT NULL 


STORE THE CORRELATION VALUE 


DO FOR THE NUMBER OF ELEMENT PLACEMENTS 


DO FOR THE NUMBER OF LINE PLACEMENTS 
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3.5.12 CORSDA (Correlation Using S.S.D.A. W Thresholding) 

This subroutine Is exactly the same as CORG with the exception 
that threshold checking is ^>er£oriDed. l/hen a landmark placement 
Is being computed, If the partial correlation value exceeds the 
previous minimum value, the correlation for this placement Is ter- 
minated and the next placement is tested. 


CORSDA 



SET THE BEST FIT ELEMENT AND LINE TO 1 


DETERMINE THE SEARCH AREA MEAN 


SET THE CLOUD COUNTER AND CORR. VALUE TO 0 


INCREMENT 


IS THIS A CLOUD? 


ADD INTO CORRELATION 


CLOUT COUNTER VALUF FOR THIS PLACEMENT 


V DOES CORRELATION VALUE EXCEED / 
\ PREVIOUS MINIMUM? A 


EXIT OUT OF THIS LOOP 


DO FOR THE // ELEMENTS IN THE LANDMARK 


DO FOR THE // LINES IN THE LANDMARK 


ADD IN CLOUD PIXELS 


IS THIS A BEST FIT? 


SAVE BEST BIT 



STORE CORRELATION VALUE 


DO FOR THE # OF ELEMENT PLACEMENTS 


DO FOR THE if OF LINE PLACEMENTS 
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3.5.13 CORSMN (Correlation Surface Main Driver) 


CORSMN Is tlie driver program to display a correlation surface. 
This routine calls three routines; GRID, SCALE, and SPLPRO. 


FORMAT: CALL CORSMN 


CORSMN 


CALL THE GRID ROUTINE 


CALL THE SCALING ROUTINE 


CALL THE DISPLAY ROUTINE 
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3.5.14 CRALGS (Correlation Algorithm Controller) 


CRALGS Is the driver to the various correlation algorithms. 

This subroutine Is called by CORALA. This subroutine will call the 
appropriate algorithm based upon the variable ITALG. ILMSZE Is the 
current landmark size and IES,ILS are the screen coordinates of the 
upper left hand comer of the current landmark. IES3 Is the search 
area's starting screen coordinate. 

Upon return from the various algorithm routines, this sub- 
routine will draw a box around the best fit placement and will store 
the actual correlation values into the array lACT. lACT will be 
stored into the data base If the data base mode Is on. 

FORMAT; CALL CRAI,GS(ILMSZE,IES,ILS .ITALG, IES3) 

CRALGS 

READ IN THE SEARCH AREA DATA 
IN AN ARRAY 

READ IN THE LANDMARK DATA 
CALL THE ROUTINE TO TURN ON CLOUD BITS 
GET THE LANDMARK MEAN 
CALL THE CORRECT CORRELATION ALGORITHM 
DRAW A BOX .GROUND THE BEST FIT 
STORE THE CORRELATION VALUES 


I 
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3.5.15 CRCORR (Create Correlation Coefficient) 

The CRCORR subroutine creates the correlation coefficient of 
the auto correlation surface versus the arcual computed correlation 
surface. This value is between 1. A good correlation is one with 
a high positive correlation coefficient. RCORR is the value of the 
correlation coefficient. Algorithm used: 


81 ^ 

I (lAUTOd) - IAUTO(I))*[lACT(D - IACT(I)J 
I»T 


81 81 \h 

I (lAUTOd) - IAUT0(I))2*[ (lACTd) - lACTd))^) 


FORMAT: CALL CRCORR (RCORR) 
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3.5.16 DISPLE (Display the Slope) 

When the user specifies that a correlation surface should be 
displayed, DISPLE will display the best fit's correlation value 
along with a one placement buffer on all sides. This displays a 
3x3 matrix of correlation values. Both the actual and the auto 
correlation values are displayed in the RM9100. 

FORMAT: C.'LL DISPLE 
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3.5.17 DISPLY (Display Imagery) 


The DISPLY subroutine displays an M x N matrix of picture ele- 
ments on the RM9100 beginning at line '>'LS and element lES. The 
array lARRAY contains the data to be displayed. This array contains 
data in an lARRAY (element //, line if) format. lELEMS specifies the 
number of elements and ILINES the number of lines to be displayed. 

FORMAT: CALL DI SPLY( lARRAY ,IBUF .lELEMS , ILINES ,IES ,ILS) 

DISPLY 


OET A LINE OF DATA IN 
lARRAY STORING INTO IBUF 

WP.ITE IBUF TO THE RM9100 

DO FOR THE NUMBER OF LINES 
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3.5.18 DIVID (Classify the Image) 

DIVID Is used to create a classi ^i.ed Image using an Image's 
.65 micron and .85 micron data. This routine will classify data to 
four categories: 1) water; 2) vegetation; 3) bare earth; and 4) 
clouds, snow, or Ice. 

A ratlolng algorithm Is used for this classification. The 
VR/IR ratio defines the classification. 


VR 


This routine will classify a square ISIZE x ISIZE area. The 
.65 and .85 micron data must currently be displayed on the RM9100. 

The point (lESl.ILS) defines the upper left comer of the .65 micron 
data and the point (IES2,ILS) the corresponding .85 micron data 
location. The point (IES3, ILS) is the upper left hand comer of 
the display area the classified image will be written to. The 
classified image is written out in A grey levels — 25, 50, 75, and 
100 corresponding to water; vegetation; bare earth; and clouds, snow, 
or ice, respectively. IFLAG flags whether the classified image is 
to be written to the RM9100. If IFLAG equals one, the data is not 
written; otherwise it is. ICLOUB is the percentage of cloud pixels 
in this scene. 

FORMAT: CALL DIVID(ISIZE, ILS ,IES1 ,IES2 ,IES3 , IFLAG ,ICL0UD) 
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DIVID 


INITIALIZE CLOUD COUt-TER 


GET THE NEXT ELEMENTS VR AND 
IR COMPONENT 


IS THE RATIO LESS THAN 
THE VEGETATION THRESHOLD? 




IS THE VR LESS THAN THE 
CLOUD THRESHOLD 



\ IS THE RATIO LESS / 

aASSIFY AS 

\ THAN THE BARE / 

VEGETATION 

EARTH THRESHOLD /p 


CLASSIFY AS 

CLASSIFY 


AS BARE 

AS 


EARTH 

WATER 


INCREMENT 

CLOUD 

COUNTER 


CLASSIFY 
AS A CLOUD 


DO FOR THE NUMBER OF ELEMENTS 


IF IFLAG 4 1 WRITE OUT THE CLASSIFIED LINE 


DO FOR THE NUMBER OF LINES 


DETERMINE CLOUD PERCENTAGE 
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3.5.19 GRID (Correlation Surface Grid) 


GRID is used with the correlation surface routines. This sub- 
routine draws on the RM9I00 a 3-D type grid on which the correlation 
surface will be plotted. The final output of GRID is the following: 



AUTO ACTUAL 


FORMAT: CALL GRID 


GRID 


RESET THE RM9100 


DRAW THE HASH MARKS 


DRAW THE X AXIS 


DRAW THE Y AXIS 


DRAW THE Z AXIS 


LABEL THE X AXIS 


I 
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3.5.20 ^.SAHIS 


This algorithm attempts to equalize the output data Intensities 
so that the entire range of grey scales are used. The range 0-127 
or 0 - 63 Is used by the output data. This program creates the 
cumulative probabllltv function and uses this for display purposes. 

FORMAT: CALL GSAHIS (ILS.IFLAG2) 

GSAHIS 


DECID'- WHICH BAND TO USE 


^ IS THIS 

BAND 1? ^ 

EQUALIZE TO 

EQUALIZE TO 

128 VALUES 

64 VALUES 


PERFORM THE EQUALIZATION 


OUTPUT THE EQUALIZED BAND 
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3.5.21 GSASTD 


This algorithm uses a standard deviation method for grey scale 
adjustment of Imagery. The data Is adjusted to 128 or 64 grey 
scale values. The data Is scaled using the standard deviation of 
the data and the number of output Intensities the user desires. 


FORMAT: CALL GSASTD (ILS,IFLAG2) 


GSASTD 


DECIDE WHICH BAND TO USE 

IS THIS BAND 1? 

EQUAI.IZE TO 

EQUALIZE TO 

128 V.UUES 

64 VALUES 

PERFORM THE 

EQUALIZATION 

OUTPUT THE EQUALIZED BAND 
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3.5.22 SAREA (Search Area Extraction) 


Thla routine allows the user to move a cursor over the scene 
displayed on the monitor and select a search area to be used for 
correlation purposes. When this Is done, the original scene Is 
erased and two representations of the area are displayed on the 
upper portion of the screen. These representations are the two 
spectral scenes. 

ILS is the starting line that the search area should be dis- 
played at once they are defined by the user. IX and lY upon Input 
are the starting coordinates for the cursor. On output, they are 
the upper left hand comer coordinates of the defined search area. 

FORMAT: CALL SAREA (ILS.IX.IY) 


SAREA 


\ IS THIS A DATA BASE STORAGE RUN? / 

T\ / F 

SEARCH AREA SIZE IS 64 

GET THE SEARCH AREA SIZE 
FROM THE USER 

IS THIS A BATCH CONTROLLED RUN? 

READ IN SEARCH AREA START 
LINE, ELEMENT 

1 

GET SEARCH AREA START LINE, 
ELEMENT FROM THE USER 

— 1 


REDRAW SCREEN WITH THE SEARCH AREAS 
BEING DISPLAYED 
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3.5 23 SCALE (Correlation Surface Scaling) 

The SCALE subroutine acalea the correlation aurface valuea to 
fit to a grid ayatem defined by the GRID aubroutlne. The GRID 
routine defines a 180 x 180 pixel x,y grid. A perspective view la 
given by defining a z axis. 



IXREF Is an array that will contain the scaled x points, and 
lYREF will contain the scaled y points of a 9 x 9 correlation surface. 
The best fit's correlation value will always be at 

[i:<REF(5,5)IYREF(5,3)) . 

The correlation surface will contain the best fit and a 4-placement 
buffer on all sides. 

SCALE 


FIND THE MAXIMUM VALUE 
TO SCALE TO 


CREATE THE X SCALE ARRAY 


CREATE THE Y SCALE ARRAY 
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3.5.24 SPIX 

This subroutine performs tenth-of-a-plxel registration of 
imagery. This routine used a cubic convolution process to perform 
the interpolation. 

The array iTWTN contains the best fit correlation value along 
with a 2>placement buffer on all aides. This array is expanded to 
a 50 X 50 element array using the convolution algorithm. The beat 
tenth-of-a-plxel fir can then be deteriined. 

FORMAT: CALL SPIX ( ITALG , ILMSZE) 

SPIX 


CREATE THE 5 X 5 ELEMENT CORRELATION SURFACE 


CREATE THE 50 X 50 ELEMENT CORRELATION SURFACE 




IS THIS THE CLASSICAL CORRELATOR ALGORITHM? 


MAXIMUM VALUE IN BEST 
SUBPIXEL FIT 


MINIMUM VALUE IN BEST 
SUBPIXEL FIT 
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3.5.25 SPLPRO (Spline Processing) 

SPLPRO Is the main processing routine for the correlation 
surface display. This routine basically takes a 9 x 9 pixel 
correlation surface and Interpolates this to a 180 x 180 pixel 
surface. The interpolation uses a cubic polynomial spline fit. 

This routine basically fills in the holes in the correlation surface. 
EXAMPLE: 

• • 

4 pt. line , , 

20 pt. interpolated » , , • ‘ . 

line * 

This routine fiist processes the horizontal splines and the 

vortical splines. 

FORMAT: CALL SPLPRO (IXREF , lYREF) 

SPLPRO 

PROCESS HORIZONTAL SPLINES 
PROCESS VERTICAL SPLINES 
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3.5.26 WEICFTr (Determine Weighting Factor) 


The WEIGHT function determines the weighting factor to be used 
in a cubic convolution subpixel registration. The algorithm used 
is as follows: 

weight • I - 21x1^ + |xP if Oi. jx] u. I 
weight ■ 4 - 8lx| + 5|xp - [x|^ if I lx| i.2 
RDEL is the x value in the equations. 

FORMAT: y - WEIGHT (RDEL) 

WEIGHT 


OKT ABSOLL^TE 

VALUE OF RDEL 

\ IS RDEL >1? y/ 

WEIGHT “ 4-8*RDEL+ 

WEIGHT - 1-2*RDEL+ 

5**RDEL**2-RDEL**3 

RDEL**3 


h(x) 
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3.6 Library Routines 


The library routines are those subroutines which perform a general 
function without being strictly I/O routines. These routines are; 

SECTION 


1. 

DECDE 

3.6.1 

2. 

DRABOX 

3.6.2 

3. 

GETHDR 

3.6.3 

4. 

IMGTT 

3.6.4 

5. 

LMDISP 

3.6.5 

6. 

MEAN 

3.6.6 

7. 

NXDREC 

3.6.7 

8. 

POS 

3.6.8 

9. 

REDAT 

3.6.9 

10. 

REDRAW 

3.6.10 

11. 

SPCOEF/SPLINE 

3.6.11 
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3.6.1 DECDE (Decode) 


The DECDE subroutine works much like a FORTRAN DECDE statement. 
This routine allows the user to decode an integer into an ASCII 
array . 

ABUF is the byte array to contain the decoded ASCII string. 

INUM is the integer to be decoded. ICNT is the maximum number of 
characters in INUM. 


FORMAT: CALL DECDE (ABUF, INUM, ICNT) 


DECDE 


1 

STRIP OFF THE RIGHT-MOST DIGIT 
AND ADD "60 TO IT TO MAKE 
IT AN ASCII NUMBER 

1 

SHIFT THE NUMBER OVER 
DIGIT TO THE RIGHT 

1 


DO ICNT TIMES 

STRIP OFF LEADING ZEROS 
MAKING THEM BLANKS 
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3.6.2 DRABOX (Draw a Box) 


The DRA.^OX routine draws a square box sized ISIZE x ISIZE. 
INTENSE defines the intensity that this box is to be written at on 
the RM9100. IX and lY define the upper left hand comer of the data 
to be enclosed by the box. IBUF is a scratch array and IBUFl is the 
array of data that the box overwrite^, , 



FORMAT: CALL DRABOX(INTENS , ISIZE , IX. lY . IBUF , IBUFl ) 

Restrictions: 0 <. ISIZE <, 62 

DRABOX 


INITIALIZE THE SCRATCH ARRAY 
TO THE OUTPUT INTENSITY 


READ THE IMAGERY TO BE Ol^RWRITTEN 
.\ND SAVE IN IBUFl 


WRITE THE BOX ON THE RM9100 
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3.6.3 GETHDR (Get Header Data) 


The GETHDR subroutine gets the header data from a LANDSAT image 
disk file. This routine requires that unit 3 has been opened as the 
unit specifying LANDSAT Imagery. This header data Is found at 
record 1. The common block HEADER contains the arrays to receive 
this header data. The following Items are retrieved: 

SID(ll) - 11 character satellite l.D. 

EXPODT(7) - 7 character exposure date 
F0RCEN(14) - 14 character format center 
SUNANG(2) ~ 2 character sun angle 
SUNANA(3) - 3 character sun azimuth 
NADIR(14) - 14 character nadir 

LBUP is a scratch array used to read the header record, 

FORMAT : CALL GETHDR CMUF) 

GETHDR 


READ 

THE 

HEADER 

RECORD 

EXTRACT 

AND 

STORE 

THE 

SID 

EXTRACT 

AND 

STORE 

THE 

EXPODT 

EXTRACT 

AND 

STORE 

THE 

FORCEN 

EXTRACT 

•\ND 

STORE 

THE 

SUNANG 

EXTRACT 

AND 

STORE 

THE 

SUN ANA 

EXTRACT 

ANE 

STORE THE 

NADIR 
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3. 6. A IMGTT (Get the Image and Display) 


IMGTT Is the driver routine to display LANDSAT imagery onto the 
RM9100. This routine allows the programmer to display an M x N pixel 
area on the RAMTEK. lES and IL5 are the starting element, line of 
the LANDSAT Imagery to be displayed. This data is stored in a 3584 
element x 2340 line data file. lESCRN and ILSCRN are the starting 
element, line screen coordinates at where this data will be written. 
NELEMS is the number of elements/lines to be written and NLINES the 
number of lines. IRATIO defines the output ratio of the displayed 
data. If IRATIO is 1, then every pixel is displayed. If IRATIO is 2. 
then every other pixel on every other line is displayed and so on. 

This gives the effect of shrinking down th 2 Imagery. IBUF and LBUF 
are scratch arrays. IFLAG is a flag specifying whether the header 
data is to be displayed on the screen. If IFLAG is greater then or 
equal to 1, this header information is not displayed. 


FORMAT : CALL IMGTT ( lES , ILS , lESCREN , ILSCRN , NELEMS .NLINES , IRATIO , IFLAG , 

IBUF, LBUF) 


IMGTT 


GET HE/VDER LATA 

DISPLAY HEADER DATA IF REQUIRED i 


GET A LINE OF IMAGE DATA 

DISPLAY A LINE OF IMAGERY 


DO FOR THE NUMBER OF LINES 
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3.6.5 LMDISP (Display Header) 


IXDISP displays the header data of a LANDSAT Image onto the ter- 
minal screen. This routine Is generally used after GETHDR, get header, 
has been called. 

All Input values required to be displayed are stored In the 
arrays contained In the common block HEADER. The variable IFLAG 
determines whether these arrays are to be printed. If IFLAG Is 
greater than or equal to 1, the arrays are not displayed. 

FORMAT: CALL UiD ISP (IFLAG) 

LMDISP 



IS FLAG >.1? ^ 


DISPLAY SID, EXPODT, FORCEN, 

NULL 

SUNANG, SUNANA, NADIR, 

HEADNG ON THE TERMINAL 


SCREEN 
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3.6.6 MEAN (Determine the Mean) 


The subroutine MEAN calculates the mean of a group of array ele- 
ments. The array lARRAY (Element, Line) contains the data values. 
lARKAY must be dimensioned as a 64 x 64 square 2-dlmenslonal array. 
RWM Is the calculated mean value. ISIZE Is the umber of elements 
and lines to be used. lES and ILS are the offsets Into lARRAY from 
where the mean Is to be taken. 

FORMAT: CALL MEAN ( lARRAY , ISIZE, RWM, ILS , IE.<=> 

MEAN 


INITIALIZE MEAN TO ZERO 


1 

1 

SUM UP NEXT ELEMENT TO 
PREVIOUS SUM 

DO 

FOR THE NUMBER OF ELEMENTS 

1 

DO FOR THE NUMBER OF LINES 


DIVIDE THE TOTAL SUM BY THE 
NUMBER OF ELEMENTS X THE 
NUMBER OF LINES 
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3.6.7 N XDREC (Next Disk Record) 

The subroutine NSDREC is used to fill up an array that will 
contain image data. The image data will be read from unit 8. This 
subroutine is generally called by IMGfT. 

Each line of imagery is composed of seven 512-byte records. The 
variable lAVAR defines what the desired line is, NRED is the number 
of records needed to extract the desired data. NELEMS is the number 
of elements desired, and lES is the starting element on the line. 

FORMAT: CALL NXDREC(IAVAR, NRED, NELEMS ,IES) 

NXDREC 


DETERMINE THE STARTING RECORD TO READ 

^ IS THE LINE OUT OF BOUNDS? y/ 

BLANK FULL 

READ A RECORD AND STORE 

DO FOR THE NUMBER OF 
RECORDS REQUIRED 

UNPACK THIS DATA INTO IBUF 
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POS (Poaitlon Cursor) 

POS is used CO allow the user Co poslclon Che cursor on Che 
RM9100 screen. This routine allows a 9 x 9 pixel cursor to be dravn 
and moved around by Che user. When switch 4 is clipped up, Chen 
down, Che user has positioned the cursor Co Che desired location, and 
control is returned to Che calling routine. 

IX and lY are the Initial locations ot the center of the cursor. 
IX and lY also receive the final values of the cursor center. ILSTRT, 
ILSTOP, lESTRT, and lESTOP define a rectangular area in which the 


cursor may move. 


RM9100 



cursor window 


FORMAT: CALL POS ( IX, lY , ILSTRT , ILSTOP , lESTRT , lESTOP) 

POS 



IS THE STICK IN DETENT? ^ 


READ THE DATA UNDER tfflERE THE CURSOR 


IS TO BE DRAWN AND SAVE 

NULL 

DRAW THE CURSOR 


REDRAW THE DATA 


DO UNTIL SWITCH 4 IS UP 
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3.6.9 REDAT (Read in Data) 


RFDAT will read a M x N area from the RM9100 Into a two- 
dimensional array lARRAY. The maximum size area is 64 x 64 pixels. 
The area is read beginning at line ILS and element lES. IBUF is a 
scratch array. lELEMS is the number of elements/ line, and ILINES is 
the number of lines to read in. 

FORMAT: CALL REDAT ( lARRAY , IBUF , lELEMS , ILINES , lES , ILS) 

REDAT 



READ A LINK OF DATA 
INTO THE SCRATCH ARRAY 


STORE THE SCRATCH ARRAY 
INTO THE FINAl. BUFFER 

DO 

FOR THE NUMBER OF LINES 
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3.6.10 REDRAW (Redraw cha Box) 


REDRAW la used in conjucclon with DRABOX. DRABOX draws a 
box around an area storing the data being overwritten by the box, 
REDRAW redraws this overwritten data. 

ISIZE Is the size of the box to be redrawn. IX and lY are the 
screen locations of the upper left hand comer of the area enclosed 
by the box. IBUF is the array containing the prestored data. 

FORMAT: CALL REDR/.W(ISI2E,IX,IY,IBirF) 

REDRAW 


DETERMINE BOX SIZE 


REDRAW 4 SIDES OF 
THE BOX 







3.6.11 SPCOEF/ 

SPLINE 

The subroutine SPCOFF end the function SPLINE are designed to 
evaluate the natural cubic Interpolatory spline fit to a set of data 
specified by the array of N nodes, XN, with corresponding function 
values In the array FN. The nodes of XN must be distinct. The spline 
Is determined in SPCOEF and is evaluated In SPLINE. SPCOEF arranges 
the nodes in Increasing order and stores this order in the array INDEX. 
Then XN(INDEX(l)) Is the smallest node, XN(INDEX(2)) the next largest, 
etc. The array XN itself Is not altered. SPCOEF then evalutates the 
array S of the N second derivatives needed to define the spline. The 
call arguments of the subroutine subprogram SPCOEF are: 

SPCOEF CN ,XN , FN , S , INDEX) 

The remainder ot the spline calculation is done by the function sub- 
program SPLINE, which has the call arguments: 

SPLINE (N . XN , FN , S , INDEX) 

Here N,XN,FN,S, and INDEX are as before; and X is the point at which 
we wish to evaluate the spline. Its value is returned by the function 
name SPLINE. The arrays .XN, IN, and INDEX must be dimensioned in the 
calling program. 

If X < XN(INDEX( 1) ) , the function is approximated by a straight 
line which passes through ' » point (XN(INDEX( 1) ) , FN(INDEX(l))) with 
a slope which agrees with the slope of the spline at that point. If 
X > XN(INDEX(N) ) , a straight line approximation is again used; this 
time it passes through (XN(INDEX(N) ) , FN (INDEX(N) ) ) and has the slope 
the spline at that point. 

If we want to approximate a function / j’ by a spline at several 
points X given a set of nodes XN , the routine SPCOEF need be called 
only once, and SPLINE is called once for each point x. 
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The SPCCEF/SPLINE program was defined in the bcok : 

NUMERICAL COMPUTING; An Introduction , by Lawrence ?. Shampine 
and Richard C. Allen, Jr., 1973, by WB. Saunders. 
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3.7 


I/O Roatlnes 


The I/O routines are those routines which are used strictly for I/O 
purposes, Th'ise routines are: 


SECTION 


1. 

ADJUST 

3.7.2. 1 

* 

INIJOY 

CM 

3. 

REi\DJY 

3. 7. 2. 3 

4. 

RESET 

3. 7.1. 1 

5. 

RMCHCK 

3. 7. 1.2 

6. 

RMERSE 

3.7. 1.3 

7. 

RMRDIM 

3.7. 1.4 

8. 

RKIEXT 

3.7. 1.5 

9. 

m^TRIM 

3. 7. 1.6 

!0. 

WIO 

3.7. 1.7 
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3.7.1 RA^f^EK 9100 I/O ROUTINES 


All RM9100 I/O routines use the WTQlo and GETADR system I/O sub- 
routines. The baslv icture for a RM9100 I/O routine is the follrwing: 

(1) The command instruction Is placed in a buffer. 

(2) The address of the command buffer and the length of the >"’.ffer is 
stored in a parameter array. This is where GETADR is used. 

(3) Depending on the function, the appropriate data from the parameter 
array is written/read to the RM9100. The operation of reading/ 
writing is performed by using a WTQIO, 

Two auxiliary routines that have been used by all RM9400 I/O sub- 
routines are WIO and RMCHCK. These routines have been developed so that 
all I/O functions will utilize them and help conserve task space. 
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3.7.1. I RESET (Reset the RM9100) 


The RESET instruction is an RM9100 system clear function 
identical to the power-on (or hard reset) initialization 
sequence. All pending operations are discarded. The refresh 
memory will be erased. 

FORMAT CALL RESET 
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3.7 1.2 RMCHCK (Check RM9100 I/O) 


The RMCHCK subroucliie checks to see If the previous RM9400 
I/O was correctly handled by the RM9100. The status words are 
contained In the common block 10. The directive status word, 
IDS, Is equal to 1 If the directive was sent successfully. The 
I/O status word, lOST(l), Is 1 If the issued I/O executed 
correctly. If either of these variables do not equal 1, the I/O 
was unsuccessful. 

FORMAT: CALL RMCHCK 
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3.7. 1.3 RMERSE (RM9100 Erase) 


The RMERSE allows the progranmer to erase a rectangular 
region on the RM9100. This region is defined by the variables 
ILSTRT, ILSTOP, lESTRT, and lESTOP. This is a software erase 
which uses the RMWRIM subroutine to write blank data into memory. 

FORMAT: CALL RMERSE (ILSTRT , ILSTOP .lESTRT, lESTOP) 
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3.7. 1.4 RMRDIM (RM9100 Read Image) 


The RMRDIM subroutine allows Che user Co read a rectangular 
area from the RM9100 memory. Up to 2000 pixels may be read in 
any one call to RMRDIM. 

ISCAN determines the primary and secondary scan directions. 
ISCAN “ 0 primary - left to right 

secondary - top to bottom 
ISCAN ■ 4 primary - top to bottom 

secondary - left to right 

IXS and lYS are the element and line values in memory at 
which the read is to begin. IDATA is the buffer in which data 
is to be filled, and ILEN is the number of pixels to be read. 

FORMAT: CALL RMRDIM( ISCAN, IXS ,IYS , IDATA, ILEN) 
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3.7. 1.5 RMTEXT (Write Text on the RM9100) 

The routine WfTEXT writes character text on the RM9100. 

In the COR.\LA program, the following variables should always 
have these values: 

lADD - 0 
IPOLAR - 0 
ICHAN - -1 
ISCALE - 0 

ISCAN is the variable defining the scan direction in which 
the text will be written. 

ISCAN ■ 0 primary - left to right 

secondary - top to bottom 
ISCAN ^ 4 primary - top to bottom 

secondary - left to right 

IXS and lYS are the upper left hand start point for the text. 

All text is written in 7 x 9 pixel fonts. The variable ICHAR is 
the address of the first character to be written. INCHAR is the 
length of the character string. INCHAR should always be an even 
number. 

FORMAT: CALL RMTEXT (lADD , IPOLAR , ICHAN , ISCAN , ISCALE , IXS , ICHAR , 

INCHAR) 
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3. 7. 1.6 RMWRIM (RM9100 Write Image) 


The RMWRIM writes to any rectangular area In the RM9400 
In word mode. Up to 2000 pixels may be written In any one call 
to RMWRIM, 

In the CORALA programs, ICHNL always equals *1, ISCALE 
always equals 0, and IFCTN always equals 0. ISCAN Is the scan 
direction to be used for the write. 

ISCAN - 0 primary - left to right 

secondary - top to bottom 
ISCAN ■ 4 primary - top to bottom 

secondary - left to right 

IXS and lYS are the element, line start where the write to 
memory is to occur. IDATA is the buffer address of the element (s) 
to be written, and ILEN is the number of pixels to be written. 

FORMAT: CALL RMWRIM( ICHNL , ISCAN, ISCALE , IFCTN , IXS ,IYS , IDATA, ILEN) 
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3. 7. 1.7 WIO (Issue Wait for Queue I/O) 


The WIO subroutine issues a POP 11/70 WTQIO system directive. 
There are two arguments passed to this routine, IFCTN and IDEV. 
IFCTN specifies the I/O function to be executed, and IDEV is the 
device in which the I/O function is to occur. 

This subroutine is used for all RM9400 I/O. The common 
block 10 contains the other parameters required by the WTQIO 
instruction. 

FORMAT; CALL WIO (IFCTN , IDEV) 

FUNCTION CODES: "410 - Write Logical block (RM9400) 

"1010 - Read Logical block (RM9400) 



3.7.2 JOYSTICK I/O ROUTINES 


The Joystick device on the POP 11/70 is terminal number TTA 
The Joystick routines read this terminal and interpret the data. 
For data to be sent from TT4:, one of the switches marked 1 
through 4 must be switched up. 
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3. 7.2.1 ADJUST (Adjust. Data) 


The ADJUST subroutine Is written In DEC MACRO assembly 
language. This routine adjusts the bits In the output data 
from the Joystick. This routine Is device dependent and called 
only by the subroutine READJY. The source code for ADJUST is 
found in the file ADJDAT.MAC. 
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INIJOY (Inltlalljsa the Joystick) 


The INIJOY subroutine assigns a logical unit nujiber to the 
joystick. This subroutine Is called once In CORALA. 

The parameter ILUN is the logical unit number to assign the 
joystick to. This unit number will be used by the routine 
READJY, to read the joystick device. 

FORMAT: CALL INIJOY (ILUN) 


B-76 



3.'/. 2. 3 READJY (Read the Joystick) 


The READJY subroutine allows the user to read the joystick. 
Four words will be returned on each call. These words will be 
the switch word. X. Y. and Z words. The X, Y, and Z words have 
a range between 0 and 1023. When the stick Is In detent, f’" 

X and Y words are 512. When the stick is up in the uppe- I"'' 
corner, the X and Y words are 0, The Z word Is set by ' i 
stick knob. The switch word is bit oriented and follows . 
convention : 

Bit Location 15 K 1 3 12 1 1 10 9 8 7 6 5 4 3 2 1 0 
Switch Location X X X X 

SWITCH WORD 

If bit 0 is set, switch 1 is up. If bit 1 is set, switch 2 
is up, and so on. There are 4 switches. 

The IWRDCT parameter must always be 4. ABUFF is the buffer 
to receive the 4 words of output. ILUN is the logical tmit 
number of the joystick as initialized by INIJOY. 

FORMAT ; CALL RE ADJY (ABUFF , IWRDCT , ILUN ) 
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4.0 DATA BASE DEFINITION 


A correlation analysis data base has been defined. This data base 
allows the user to store correlation run data. When the user specifies 
that the data base will be updated, certain parameters are selected for 
the user automatically. The test case number Is automatically selected 
as well as the search area size of 64; landmark sizes of 48, 32, 24, 16, 
and 8; and the selection of subpixel registration on these landmarks. 

The test Is run In band 1, .65 micron, and band 2, .85 micron. 

The data base file is located in [H2,5] CORRELATE.DAT. This Is a 
direct access file, 3584 bytes per record. The first record is a header 
record. Element 1 of the header recoro contains the current number of 
tests that have been run. Subsequent records contain test data. Each 
fest creates two records of data. One record for the test run in band 1 
and one record for the tesc run in band 2. 

The test data recorded is scored ac integers or real values dependent 
upon their location in the record. The record definition is as follows: 

CORRELATION ANALYSIS DATA BASE DEFINITION 
Two records are written per test case. There is one header record 
located at record 1, records 2 and 3 contain data for test 1, records 4 
and 5 contain data for test 2, and so on. . . 
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Correlation Analysis Data Base Definition (Continued) 


Description // of 

Words 

LM Size 

Start Buffer 

Addr . 

TEST CASE NUMBER 

1 


1 

CORFIN 

S.ID 1 2 CHARS /WORD 

6 


2 

COPINI 

S.ID 2 2 CHARS /WORD 

6 


8 


SCENE 1 START LINE 

1 


14 

CORINI 

SCENE 1 START ELEMENT 

1 


15 


SEARCH AREA 1 START LINE 

1 


16 

CORINI 

SEARCH AREA 1 START ELEM 

1 


17 


S.ID 1 LM % CLOUD *10 

1 

8 

18 

CORFIN 


1 

16 

19 



1 

24 

20 



1 

32 

21 



1 

48 

22 


S.ID 2 S.A. % CLOUD*10 

1 

1 


23 

CORINI 

S.ID 2 LM % CLOUD *10 

1 

8 

24 

CORFIN 


1 

16 

25 



1 

24 

26 



1 

32 

27 



1 

48 

28 


CORRELATION TIME SSDA*10 

1 

8 

29 

CORFIN 


1 

16 

30 



1 

24 

31 



1 

32 

32 



1 

48 

33 


CORRELATION TIME CLASS* 10 

1 

8 

34 

CORFIN 


1 

16 

35 



1 

24 

36 



1 

32 

37 



1 

48 

38 


SSDA ACT VS. AUTO Qf" . 

1 

8 

39 

CORFIN 

C0EFFICIENTS*10000. 

1 

16 

40 



1 

24 

41 



1 

32 

42 



1 

48 

43 


CLASSICAL ACT. VS. AUTO 

] 

8 

44 

CORFIN 

CORR. COEFFICIENTS* 10000. 

1 

16 

45 



1 

24 

46 



1 

32 

47 



1 

48 

48 


L.M. ST/J^T LINE S . ID 2 

1 


49 

CORFIN 
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5 . 0 AUXILIARY PROGRAMS 


Auxiliary p ograms were written to aid the user to control the flow 
of the CORALA program. These programs are: 

(1) INICORALa - Initialize the CORALA data base 

(2) AUTOCOR - This program creates the batch mode run stream for CORALA 

(3) LANDSAT - This program reformats LANDSAT tapes to a deinLer leaved 

format 

(4) FULL - Displays a 1024 x 1024 section from a reformat LANDSAT image 
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